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ABSTRACT 

We analyze archival Lya spectra of 12 “Green Pea” galaxies observed -with the Hubble Space 
Telescope, model their Lya profiles ■with radiative transfer models, and explore the dependence of 
Lya escape fraction on various properties. Green Pea galaxies are nearby compact starburst galaxies 
with [OIIIJA5007 equivalent widths of hundreds of A. All 12 Green Pea galaxies in our sample show 
Lya lines in emission, with a Lya equivalent width distribution similar to high redshift Lya emitters. 
Combining the optical and UV spectra of Green Pea galaxies, we estimate their Lya escape fractions 
and find correlations between Lya escape fraction and kinematic features of Lya profiles. The escape 
fraction of Lya in these galaxies ranges from 1.4% to 67%. We also find that the Lya escape fraction 
depends strongly on metallicity and moderately on dust extinction. We compare their high-quality 
Lya profiles with single HI shell radiative transfer models and find that the Lya escape fraction anti¬ 
correlates with the derived HI column densities. Single shell models fit most Lya profiles well, but 
not the ones with highest escape fractions of Lya. Our results suggest that low HI column density 
and low metallicity are essential for Lya escape, and make a galaxy a Lya emitter. 


1. INTRODUCTION 

The Lya emission line is a key tool in discovering and 
studying high redshift galaxies, and a good probe for 
reionization. High redshift Lya emission line galaxies 
(LAE) have been found routinely for almost two decades 
(e.g. Dey et al. 1998; Hu et al. 1998; Rhoads et al. 
2000; Ouchi et al. 2003; Gawiser et al. 2006; Guaita 
et al. 2010; Clement et al. 2012; Shibuya et al. 2012; 
Matthee et al. 2014). These high redshift LAEs are 
generally small, young star forming galaxies. They have 
compact size, low stellar mass, low dust extinction, low 
metallicity, and high specific star formation rate (sSFR) 
(e.g. Malhotra 2012; Bond et al. 2010; Gawiser et al. 
2007; Pirzkal et al. 2007; Finkelstein et al. 2008; Pen- 
tericci et al. 2009). At 2 < z < 6, these LAEs are an 
important population of star-forming galaxies, and they 
constitute an increasing fraction of Lyman break galax¬ 
ies across that range, reaching ~60% of Lyman break 
galaxies (LBGs) at redshift z~6 (Stark et al. 2011). At 
z > 6, the Lya luminosity, Lya equivalent width (EW), 
and spatial clustering of LAEs and LBGs are used to 
probe reionization of the intergalactic medium (e.g. Mal¬ 
hotra & Rhoads 2004; Ouchi et al. 2010; Hu et al. 2010; 
Kashikawa et al. 2011; Ota et al. 2010; Treu et al. 2012; 
Pentericci et al. 2014; Tilvi et al. 2014). 

To use Lya as a powerful probe of high redshift galax¬ 
ies and reionization, we should ideally understand how 
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Lya escapes from galaxies (e.g. Dijkstra et al. 2014). 
In star-forming galaxies, Lya photons come from the re¬ 
combination of gas surrounding hot young O and B stars. 
Because Lya is resonantly scattered, the path length and 
dust extinction of Lya photons is increased in a manner 
that depends on the kinematics of the gas. This process 
also determines the characteristic profile of Lya lines. 
Only a fraction of the intrinsic Lya photons can escape 
the galaxy and be observed. The emergent Lya emission 
depends on the amount of dust, the HI gas column den¬ 
sity (Nhi), the velocity distribution of HI gas, and the 
geometric distribution of HI gas and dust (e.g. Neufeld 
1990; Chariot & Fall 1993; Ahn et al. 2001; Verhamme 
et al. 2006; Dijkstra et al. 2006). 

Since the scattering of Lya photons can signihcantly 
modify the Lya profile, studying Lya proHles is an im¬ 
portant way to understand Lya escape. Lya emission 
line in LAEs usually shows an asymmetric or a double- 
peaked profile (e.g. Rhoads et al 2003; Kashikawa et al. 
2011; Erb et al. 2014). For LAEs with detected opti¬ 
cal emission lines and systemic redshifts, the peaks of 
Lya profiles are usually redshifted with respect to sys¬ 
temic velocities (McLinden et al. 2011, 2014; Chonis 
et al. 2013; Hashimoto et al. 2013; Song et al. 2014; 
Shibuya et al. 2014; Erb et al. 2014). The velocity off¬ 
set of Lya emission line from systemic velocity is usually 
smaller in LAEs than in continuum selected galaxies with 
weaker Lya emission lines or Lya absorption (Shapley et 
al. 2003). 

To understand Lya profile, the single shell outflow 
model is popular for its simplicity and ability to cap¬ 
ture several essential features of Lya lines. In a galaxy 
with spherical HI gas outflow, Lya photons backscat- 
tered from the far side of the receding shell will acquire 
a frequency shift that allows them to pass through other 
gas in the galaxy, including the near side of the shell. 
Many studies simulate the Lya radiative transfer pro¬ 
cess assuming a single such shell (e.g. Ahn et al. 2001; 
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Verhamme et al. 2006; Dijkstra et al. 2006; Schaerer 
et al. 2011; Gronke et al. 2015). The output Lya pro¬ 
file depends on the HI column density of the gas shell, 
the velocity of the shell, the dust optical depth, and the 
temperature of the HI gas in the shell. The models can 
reproduce the observed Lya profile and infer properties 
such as Nhi and outflow velocity of the HI shell. In these 
models, the Lya velocity offset is formed from scattering 
of Lya photons by gas outflows, and increasing Nhi or 
outflow velocity will usually result in larger Lya velocity 
offset. 

To get a better understanding of Lya escape, we 
need to observe high quality Lya profiles, and deter¬ 
mine systemic redshifts, gas outflows, the HI gas distri¬ 
bution/kinematics, and many other galactic properties 
of LAEs. At high redshift, however, absorption by the 
intergalactic Lya forest prevents reliable measurements 
of the blue portion of Lya emission lines. Other crucial 
observations are also impractical, both because high- 2 ; 
LAEs are faint, and because some features (notably rest- 
optical emission lines) are redshifted to \obs > ‘2A^m, 
where presently available instruments lack sensitivity. 
Therefore many studies seek to solve the Lya escape 
problem in the nearby universe by observing galaxies 
with similar properties to high-z LAEs (e.g. Giavalisco 
et al. 1996; Kunth et al. 1998; Mas-Hesse et al. 2003; 
Deharveng et al. 2008; Finkelstein et al. 2009; Atek et 
al. 2009; Scarlata et al. 2009; Leitherer et al. 2011; 
Heckman et al. 2011; Gowie et al. 2011; Wofford et al. 
2013; Hayes et al. 2005, 2014; Ostlin et al. 2009, 2014; 
Pardy et al. 2014; Rivera-Thorsen et al. 2015). How¬ 
ever, these samples of nearby Lya emission line galaxies 
have on average much smaller Lya equivalent widths and 
Lya escape fractions (1%-12% in Wofford et al. (2013) 
sample; Hayes et al. 2014) than do high-z LAE samples 
(>17% Zheng et al. 2012; Nakajima et al. 2012). The 
main exception is the present Green Pea galaxy sample. 

Green Pea galaxies were discovered in the citizen sci¬ 
ence project Galaxy Zoo, in which public volunteers mor¬ 
phologically classified millions of galaxies from the Sloan 
Digital Sky Survey (SDSS). These are compact galax¬ 
ies that are unresolved in SDSS images. The green 
color is because the [OHI] doublet dominates the flux 
of SDSS r-band which is mapped to the green chan¬ 
nel in the SDSS’s false-color ^n-band images. They 
have redshifts 0.11 < 2 : < 0.36, small stellar masses 
~ 10® — IO^^Mq, low metallicities, high specific star for¬ 
mation rates (sSFR), and emission line equivalent widths 
(EW(Ha) and EW([OHI]A5007 ) exceeding hundreds of 
A (Gardamone et al. 2009; Izotov et al. 2011). Thus 
Green Peas are good counterparts to high-redshift LAEs 
in size, morphology, stellar mass, metallicities, and op¬ 
tical emission line strengths. This suggested that Green 
Peas might also show strong Lya emission lines. 

In this paper, we use archival HST Lya spectroscopy of 
Green Peas to show that Green Peas are the best analogs 
of high-z LAEs in local universe. With high quality Lya 
spectra and rest-frame optical spectra. Green Peas pro¬ 
vide a good opportunity to study Lya escape. We explore 
relations of Lya escape to Lya profiles and to galac¬ 
tic properties, compare the Lya proHles with radiative 
transfer models, and discuss constrains on HI gas and 
Lya escape. 


2. GREEN PEAS SAMPLE AND SPECTRA DATA 
2.1. Green Peas Sample 

While the full SDSS data set contains a few tens of 
thousands of Green Peas candidates based on photo¬ 
metric selection, they are not in a category systemati¬ 
cally targeted for spectroscopic followup. In SDSS DR7, 
a sample of about 251 Green Peas were observed as 
serendipitous spectroscopic targets (Gardamone et al. 
2009). Which Green Peas were thus targeted was es¬ 
sentially a random process, depending where the SDSS 
fibers were undersubscribed by other objects of all types. 
Thus, we do not expect this step of the selection to intro¬ 
duce any important bias in the sample. A subset of these 
objects have sufficient signal to noise ratio (S/N) in both 
continuum and emission lines (Ha, H/3, and [OHIJA5007 
) to study galactic properties such as SFR, stellar mass, 
and metallicity. Galaxies with an active galaxies nucleus 
(AGN) (diagnosed by their broad Balmer emission lines 
or Ha/[NH] vs. [OHI]/H/3 diagram) are excluded. These 
selections result in 66 Green Peas that have good opti¬ 
cal spectra and measured galactic properties (see Carda- 
mone et al. 2009 and Izotov et al. 2011 for details about 
selection of Green Peas sample). This provides a parent 
sample of Green Peas for Lya emission studies. 

We searched for UV spectra of the parent sample of 66 
Green Peas in the HST archive, and find 12 Green Peas 
that have UV and Lya spectra taken with the Gosmic 
Origins Spectrograph (GOS) (Pis: Henry (GO: 12928); 
Jaskot (GO: 13293); Heckman (GO: 11727)). We study 
this sample of 12 Green Peas in this paper. The 9 galax¬ 
ies in Henry’s sample were selected from the parent sam¬ 
ple of 66 Green Peas by their FUV brightness, with 
'mpuv < 20 (AB). The 2 galaxies in Jaskot’s sample 
were selected by their extreme [OHI]/[OH] ratios. The 
one galaxy in Heckman’s sample was selected by its high 
FUV luminosity, high UV flux, and compact size. Gom- 
pared to the parent sample, the current sample of 12 
Green Peas covers the full ranges of mass, SFR, and 
EW([OHI]) of the parent sample, and these 12 Green 
Peas are only slightly biased to lower metallicity and 
lower dust extinction (see figure 1). To address the bias 
of the current 12 Green Peas sample and explore the 
relations between Lya and galactic properties, we are 
taking Lya spectra of a larger sample of Green Peas to 
cover more completely the distributions of metallicities 
and dust extinctions (PI: Malhotra (GO: 14201)). 

2.2. UV and Optical Spectra 

12 Green Peas were observed by HST-COS. All were 
imaged using the GOS acquisition modes AGS/IMAGE, 
which took high resolution NUV acquisition images and 
centered the target accurately (error ~ 0.05 arc-second) 
in the 2.5 arc-seconds diameter Primary Science Aper¬ 
ture. These images show the NUV sizes of Green Peas 
are compact compared to the 2.5 arc-seconds aperture. 
Rest-frame spectral coverages are roughly ~ 950-1500A. 
Lya spectra were taken with both FUV grating G130M 
and G160M for two sources. Since the Lya were detected 
at very high S/N ratio, we only use G160M spectra for 
these two sources. 9 sources have Lya spectra taken only 
with G160M. One source has Lya spectra taken only with 
G130M. 

We retrieved COS spectra for these 12 Green Peas 
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TABLE 1 
Sample 


ID 

RA 

DEC 

Redshift 

E(B-V)mw 

Lya flux 

EW(LyQ) 

nLya 

J esc 

12-elog(0/H) 

E(B-V) 


J2000 

J2000 


mag 

10“^^ erg cm~^ 

A 



mag 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

GP1457-I-2232 

14:57:35.13 

-1-22:32:01.8 

0.148611 

0.0410 

0.46 

7.36 

0.014 

8.04 

0.061 

GP0303-0759 

03:03:21.41 

-07:59:23.2 

0.164880 

0.0845 

1.01 

7.20 

0.050 

7.86 

0.000 

GP1244-I-0216 

12:44:23.37 

-1-02:15:40.5 

0.239426 

0.0211 

1.89 

39.96 

0.065 

8.17 

0.062 

GP1054-I-5238 

10:53:30.83 

-1-52:37:52.9 

0.252638 

0.0126 

1.54 

10.66 

0.068 

8.10 

0.069 

GP1137-I-3524 

11:37:22.14 

-1-35:24:26.7 

0.194390 

0.0156 

3.81 

33.43 

0.130 

8.16 

0.043 

GP0911-I-1831 

09:11:13.34 

-1-18:31:08.2 

0.262200 

0.0243 

3.15 

49.53 

0.155 

8.00 

0.168 

GP0926-I-4428 

09:26:0.44 

-1-44:27:36.5 

0.180690 

0.0156 

6.36 

40.82 

0.245 

8.01 

0.074 

GP1424-I-4217 

14:24:05.73 

-1-42:16:46.3 

0.184788 

0.0087 

8.55 

78.27 

0.266 

8.04 

0.028 

GP0815-I-2156 

08:15:52.00 

-1-21:56:23.6 

0.140950 

0.0352 

4.01 

75.09 

0.299 

8.00 

0.014 

GP1133-I-6514 

11:33:03.80 

-1-65:13:41.3 

0.241397 

0.0093 

2.08 

35.29 

0.352 

7.97 

0.040 

GP1249-I-1234 

12:48:34.64 

-1-12:34:02.9 

0.263389 

0.0256 

5.28 

94.80 

0.383 

8.11 

0.084 

GP1219-I-1526 

12:19:03.98 

-1-15:26:08.5 

0.195599 

0.0224 

13.45 

157.54 

0.672 

7.89 

0.000 


Note. — Column Descriptions: (1) Object ID; (4) The Milky Way extinction E{B — V)mw based on Schlafly & Finkbeiner (2011); (6) Lya flux 
measured from the spectra, and corrected for Milky Way dust extinction using the Fitzpatrick (1999) extinction law; (7) rest-frame Lya equivalent 
width; (8) Lya escape fraction; (9) metallicity measured from direct Tg method by Izotov et al. 2011; (10) Dust extinction of Green Pea galaxies. 
Errors of these measurements are dominated by systematics, so statistical errors are not given. 


TABLE 2 

Lyq; Profile Analysis 


ID 

V(blue-peak) 

V(red-peak) 

V(valley) 

FWHM(Ho) 

FWHM (red) 

Flux{blue) 

Flux(red) 

fx(valley) 
fx {red—peak) 


km 

km 

km 

km 

km 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

GP1457-I-2232 

-329 

406 

65 

188 

350 

0.410 

0.013 

GP0303-0759 

-313 

153 

-106 

233 

273 

0.047 

0.023 

GP1244-I-0216 

-240 

247 

-21 

233 

268 

0.348 

0.005 

GP1054-I-5238 

-266 

192 

-1 

291 

230 

0.108 

0.015 

GP1137-I-3524 

-355 

201 

-27 

287 

253 

0.120 

0.022 

GP0911-I-1831 

-278 

81 

-62 

298 

224 

0.185 

0.033 

GP0926-I-4428 

-223 

244 

-38 

280 

307 

0.154 

0.071 

GP1424-I-4217 

-150 

224 

78 

259 

187 

0.631 

0.059 

GP0815-I-2156 

-121 

144 

-15 

193 

185 

0.432 

0.115 

GP1133-I-6514 

-69 

271 

125 

197 

219 

0.399 

0.266 

GP1249-I-1234 

_a 

83 

_a 

209 

322 

_a 

_a 

GP1219-I-1526 

-76 

176 

25 

244 

177 

0.389 

0.114 


Note. — Column Descriptions: (1) Object ID; (2) V(blue-peak) is the velocity of the peak in the blue side of Lyo: profile; (3) V(valley) is the 
velocity of the center valley of Lya profile; (4) V(red-peak) is the velocity of the peak in the red side of Lya profile; (5) FWHM of Ha emission line 
from SDSS spectra; (6) FWHM(red) is the FWHM of Lya red peak; (7) Flux(blue)/Flux(red) is the flux ratio of blue and red part of Lya profile; 
(8) fx{valley)/f\{red — peak) is the flux density ratio of central valley to red peak of Lya profile. All velocities are relative to the peak of the Ha 
emission. 

^ GP1249-hl234 doesn’t show double-peaked Lya profile. 


from the HST MAST archive after they have been pro¬ 
cessed through the standard COS pipeline CALCOS ver¬ 
sion 3.0 (2014-10-30). The resulting spectral resolu¬ 
tions are coarser than the point-source spectral resolu¬ 
tion (FWHM~20 km s“^), because these Green Peas are 
resolved in their NUV acquisition images. Their real 
spectral resolution depends on source angular sizes. For 
FUV continuum the resolutions are about 20-50 km s“^ 
based on their NUV sizes (James et al. 2014; Henry et 
al. 2015). For Lya emission line spectra, the resolution 
may be somewhat worse if Lya is more extended than 
the UV continuum emission. If Lya filled the aperture 
uniformly, the spectral resolution would be ~200 km s“^ 
FWHM (France et al. 2009). We bin the reduced COS 
spectra to ~ 0.12A pixel~^ for 11 spectra taken with 
G160M and to ~ O.OSA pixel~^ for the one spectra taken 
with G130M. The velocity precisions of the UV spectra 
are better than 40 km s“^ (Henry et al. 2015). The pre¬ 
cisions of the systemic redshifts from Ha emission lines in 
SDSS spectra are better than 20 km s“^. The resulting 


Lya line profiles are shown in figure 2. 

Then we measure the properties of Lya emission 
lines. In most cases, we estimate the continuum level 
from rest-frame wavelength range ~ 1225-1260A. For 
GP1457-I-2232 (see Table 1 for the source ID), the Lya 
spectrum shows damped absorption wings, and we mea¬ 
sure the continuum from rest-frame wavelength range ~ 
1270-1300A. For GP1137-I-3524, we don’t have spectra 
red-ward of Lya, so we measure the continuum from rest- 
frame wavelength range ~ 1180-1210A. To get the Lya 
line flux, we integrate the continuum subtracted spec¬ 
tra in rest-frame wavelength range ~ 1210-1220A. For 
GP1457-I-2232 with damped absorption, we get its resid¬ 
ual Lya emission line flux in the center of damped ab¬ 
sorption without subtracting the continuum. We correct 
for Milky Way extinction using the attenuation deter¬ 
mined by Schlafly & Finkbeiner (2011) and the Fitz¬ 
patrick (1999) extinction law. The foreground extinc¬ 
tion values were obtained from the NASA/IPAG Galactic 
Dust Reddening and Extinction tool. 
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Fig. 1. — Green Pea sample distribution in the diagram of Hq;/H/ 3 
and metallicity. Ho/H/3 is the flux ratio of Ho and H/3 emission 
lines. In the absence of dust reddening, we expect Hq;/H/3p=j 2.86, 
and higher values of the ratio indicate the presence of dust. Red 
dots are the current sample of 12 Green Peas with Lyo spectra. 
Blue empty squares are the parent sample of 66 Green Peas. 

We downloaded optical spectra of these Green Peas 
from the SDSS DR12 archive, along with the pipeline 
measurements of their emission line properties. We cor¬ 
rect the measured Ha and H/3 fluxes for Milky Way ex¬ 
tinction (again using the Fitzpatrick (1999) extinction 
law). Then we calculate E(B-V) from dust in Green 
Peas galaxies using Calzetti et al. (2000) extinction law 
and an intrinsic Ha/H/3 ratio of 2.86, and correct dust 
reddening of the observed Ha flux. We use the metal- 
licities and mass measured from SDSS spectra by Izotov 
et al. (2011). The metallicities were calculated using 
Tf. method. These properties of Green Pea galaxies and 
their Lya emission lines are shown in Table 1. 

Then we calculate the rest-frame EW(Lya), 
EW (Lya)=Flux(Lya) //a {continuum) / (l-|-redshift), 
and the Lya escape fraction is defined as 

the ratio of observed Lya flux to intrinsic Lya flux. 
Assuming case-B recombination, the intrinsic Lya flux 
is 8.7 times dust extinction corrected Ha flux. Thus the 
/isc“ is Lya/(8.7xiJacorrected)- The SDSS Ha spectra 
were taken with 3 arc-seconds diameter aperture which 
matches the COS 2.5 arc-seconds diameter aperture 
very well. 


3. GREEN PEAS ARE ANALOGS OF HIGH-Z LAES 

The first remarkable result is all 12 Green Peas show 
Lya emission lines (figure 2; and see also Henry et al. 
(2015) for 10 of these 12 Green Peas). Furthermore, 
9 of 12 Green Peas have EW(Lya) larger than 30 A, 
and would be selected as LAEs in high redshift sam¬ 
ples. This is in contrast to other low-redshift samples 
of star-forming galaxies and Lya galaxies. The LARS 
sample (Ostlin et al. 2014), selected on the basis of high 
EW(Ha), shows Lya in emission in 12 of 14 galaxies, 
but only 6/14 would be picked up as high-z Lya emit¬ 
ters on the basis of EW(Lya) alone (Hayes et al. 2014). 
Kunth et al. (1998) detected Lya emission in 4 of 8 low- 
z galaxies in their sample, with a mean EW(Lya) of 26 


A among the detections. Wofford et al. (2013) detected 
Lya emission in 7 of 20 low-z star-forming galaxies, and 
all detections have EW(Lya)<12 A. The GALEX z~0.3 
Lya emitter sample, selected with EW(Lya)>15A, has 
mean EW(Lya) of 30 A (Deharveng et al. 2008; Cowie 
et al. 2011; Finkelstein et al. 2009; Scarlata et al. 2009). 

We then compare the EW(Lya) distribution of Green 
Peas with nearby and high-z LAE samples (figure 3). 
For the high-z samples, we use a z = 2.8 narrow-band 
selected LAE sample (Zheng et al. in preparation) and a 
sample of spectroscopically confirmed LAEs at z=5.7 and 
6.5 (Kashikawa et al. 2011). All known AGNs in these 
high-z samples are excluded, and the AGN contamina¬ 
tion for high-z LAE samples is less than 5% (Zheng et 
al. 2013). The EW(Lya) of high-z samples are com¬ 
puted from deep narrow-band and broad-band photome¬ 
tries. In figure 3, we compare the cumulative normalized 
EW(Lya) distribution for these samples. The EW(Lya) 
distribution of Green Peas is very similar to the high red- 
shift (z = 2.8) sample. According to K-S test results, the 
EW(Lya) distribution of the GALEX LAE sample dif¬ 
fers from that of the z = 2.8 LAE sample at a probability 
of 99%, while the EW(Lya) distribution of Green Peas 
sample is not distinguishable from that of the z = 2.8 
LAE sample . 

Of the two Green Peas (GP1457-I-2232 and 
GP0815-I-2156) selected by extreme [OHI]/[OH] ra¬ 
tio, one shows weak Lya and the other shows strong 
Lya emission (figure 2). They are not especially dis¬ 
tinguishable in their or EW(Lya) from the other 

Green Peas. And we don’t see correlations between 
and [OHI]/[OH] ratio in this sample of 12 Green Peas. 
So including the two sources with extreme [OHI]/[OH] 
ratio doesn’t bring obvious bias to this sample. 

4. Lya PROFILES ANALYSIS & Lyo ESCAPE 

We now use Green Peas to study the mechanisms of 
Lya escape. The escape of the resonantly scattered Lya 
photons involves potentially complex interactions with 
gas and dust. This results in richly varied Lya line pro¬ 
files that carry a great deal of information, especially 
when the profiles are double-peaked as they tend to be 
for the Green Peas. We therefore begin by exploring re¬ 
lations between Lya escape fraction and Lya profiles. 

In figure 2, we display the Lya profiles of these 12 
Green Peas in order of increasing Lya escape fraction 
from left to right, and top to bottom. The zero point of 
the velocity scale is set by the peak of the Ha line. We de¬ 
fine the “red peak” as the peak in the Lya line profile oc¬ 
curring at velocity >0. 11 of the 12 galaxies have double- 
peaked line profiles, and for these we define the “blue 
peak” as the Lya peak at velocity < 0, and the “valley” 
as the flux minimum between the two peaks for double 
peak profiles. Simply looking at their profiles, we can 
find that as increases, the residual flux at the inter¬ 
peak valley gets stronger, and the blue peak moves nearer 
to the systemic velocity. To characterize those trends 
quantitatively, we measure the velocity and flux density 
fx at the red peak, blue peak, and valley. We also mea¬ 
sure the integrated flux and full width at half maximum 
(FWHM) for the red (V>0) and blue (V<0) portions 
of the Lya profile - Flux(red), Flux(blue), FWHM (red), 
and FWHM(blue). Then we investigate correlations be- 
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Fig. 2.— Lyo and Ho emission line profiles of Green Peas. Lyo and Ho are in the same flux units of 10 erg cm ^ s ^ {km s ^ 

for all plots. These 12 galaxies are sorted by increasing from left to right, and top to bottom. The is given in each panel. 


tween those quantities and Lya escape fraction. 

As shown in figure 4, strongly correlates with 

the blue peak velocity V(blue-peak) (see also Henry et 
al. 2015) for the 11 Green Peas with double peaked pro¬ 
files. The Spearman rank coefficient for this correlation 
is r=0.84 with a a null probability of P=0.001 (Table 2). 

As the HI gas that absorbs/scatters Lya may have 
a velocity distribution similar to the HII regions in the 
galaxy, and increasing the velocity dispersion of the HI 
gas may result in larger velocity offset (e.g. Steidel et 
al. 2010), we normalize the blue peak velocity by the 
FWHM of Ha. The shows an even tighter cor¬ 

relation with the ratio V(blue-peak)/FWHM(Ha), with 
Spearman rank coefficient r=0.94 and null probability 
P=le-5. Thus adding a new variable FWHM(Ha) im¬ 
proves the correlation. As the FWHM(Ha) of these 
Green Peas have a small range of about 180-300 km s“^, 


we do a simple simulation to test the significance of the 
improvement. In each run we substitute a random num¬ 
ber between 180-300 km s“^for the FWHM(Ha) of these 
Green Peas, and calculate the Spearman rank correla¬ 
tion between and V(blue-peak)/FWHM(Ha). We 
run it 10^ times, and 1.8% of the resulting correlation 
coefficients are better than the observed data (r=0.94). 
So there is 98.2% probability that the improvement of 
correlation is not a random event. 

The doesn’t correlate with the red peak ve¬ 

locity (Spearman r=-0.2, P=0.55) or the ratio V(red- 
peak)/FWHM(Ha) (Spearman r=0.04, P=0.9), while 
it does anti-correlate with the peak separation V(red- 
peak)-V(blue-peak) (Spearman r=-0.83 and P=0.002) 
and (V(red-peak)-V(blue-peak))/FWHM(Ha) (Spear¬ 
man r=-0.75 and P=0.007). Thus correlations of 
feac^ with velocity differences V(red-peak)-V(blue-peak) 
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EWdya) [A] 


Fig. 3.— Here we compare the rest-frame EW(Lya) distribution 
of Green Peas with different samples. The solid green line shows 
the 12 Green Peas. The blue dash-dot line shows the GALEX 
z=0.3 LAE sample (Cowie et al. 2011; Finkelstein et al. 2009; 
Scarlata et al. 2009). The dashed magenta line shows the z=2.8 
LAE sample from Zheng et al. (2016, in preparation). The dotted 
red line shows the z=5.7 and 6.5 LAE sample from Kashikawa et 
al. (2011). The EW(LyQ:) distribution of Green Peas is similar to 
the z=2.8 LAE sample. 

are mostly due to the stronger correlation with blue 
peak velocity. We also notice that in figure 2 Green 
Peas with large (GP1424+4217, GP1133+6514, 

GP1219+1526) show V(Valley)>0, while the others with 
lower have V(Valley)<0. 

In figure 5, the shows a very weak correlation with 
the fiux ratio of blue and red part, Flux(blue)/Flux(red) 
(Spearman r=0.44, P=0.17). As the increases, both 
the red and blue peaks get stronger relative to the intrin¬ 
sic Lya flux. There is only a weak trend that the blue 
peak increases faster with than the red peak. In the 
middle panel of figure 5, the shows correlation with 
the residual flux density f\{valley) normalized by /a( red 
peak) (Spearman r=0.89 and P=2e-4). In the right panel 
of figure 5, shows anti-correlation with FWHM(red) 

(Spearman r=-0.83 and P=0.002). It is not clear whether 
fesi^ also has weak correlation with FWHM(blue), as 
the blue part profiles are noisy and the uncertainties of 
measured FWHM(blue) are too large. The doesn’t 
correlate with FWHM(Hq;) (Spearman r=0.04, P=0.9). 

Although the current sample size is small, these corre¬ 
lations are encouraging. We interpret these correlations 
between Lya escape and Lya profile in section 7.1. 

5. Lya ESCAPE, DUST AND METALLICITIES 

In order to probe what factors are dominant in allowing 
Lya to escape, we also investigate relations between 
and galactic properties measured from optical spectra, 
namely the dust reddening E(B-V) and the metallicity. 
In the Lya escape process, dust can absorb the Lya pho¬ 
tons and decrease f^sc^- Multiple studies have confirmed 
that Lya escapes more easily from galaxies with low dust 
reddening, and there is an anti-correlation between 
and E(B-V), although dust extinction is not the only fac¬ 
tor to determine (Atek et al. 2008, 2014; Scarlata 
et al. 2009; Finkelstein et al. 2011; Cowie et al. 2011; 


TABLE 3 

Relations of /iic“ and other, properties 




Variables 

Spearman r 

P 



(1) 

(2) 

(3) 

cLya. 
J esc 

VS. 

FWHM(Ha) 

0.04 

0.9 

pLya. 
J esc 

VS. 

V(red-peak)/FWHM(Ha) 

0.04 

0.9 

pLyce 
j esc 

vs. 

V(red-peak) 

-0.2 

0.55 

pLya. 
J esc 

vs. 

Flux(blue)/Flux(red) 

0.44 

0.17 

pLyct 
J esc 
pLya 
J esc 

vs. 

vs. 

E(B-V) 

V {red —peak) —V {blue —peak) 

-0.62 

-0.75 

0.05 

0.007 

FWHM{Ha) 

pLya 
J esc 

vs. 

metallicity 

-0.82 

0.004 

pLya 
J esc 

vs. 

FWHM(red) 

-0.83 

0.002 

pLya 
J esc 

vs. 

V (red-peak)-V (blue-peak) 

-0.83 

0.002 

pLya 
J esc 

vs. 

V(blue-peak) 

0.84 

0.001 

pLya 
J esc 

vs. 

/A(valley)//A(red peak) 

0.89 

0.0002 

pLya 
J esc 

vs. 

V(blue-peak)/FWHM(Ha) 

0.94 

10“® 


Note. — Column Descriptions: (1) Relations between and 

other properties; (2) The Spearman rank coefficient for each correla¬ 
tion; (3) The null probability of each correlation. 


Hayes et al. 2014). In the left panel of figure 6, we show 
the - E(B-V) relation for Green Peas. It is qual¬ 
itatively similar to previous studies, though this Green 
Peas sample has higher and only covers a small 

E(B-V) range. There are two outliers in figure 6. One 
is GP0303-0759 (marked by a large circle) with the low¬ 
est metallicity and E(B-V). Its WISE three bands colors 
suggest it an AGN candidate (Malhotra et al. in prepara¬ 
tion). The other one is GP1249-I-1234 (marked by a large 
square) which has single peaked line profile and large 
E(B-V). After both GP0303-0759 and GP1249-L1234 are 
excluded, the anti-correlation between and E(B-V) 
has Spearman coefficient r=-0.62 and P=0.05 (note that 
r=-0.35 and P=0.28 if only the AGN candidate GP0303- 
0759 is excluded). It suggests dust extinction is an im¬ 
portant factor of Lya escape even in this Green Peas 
sample covering a small E(B-V) range. 

Dust reddening is usually higher in galaxies with higher 
metallicity and stellar mass. So we expect f^y°‘ may 
also anti-correlate with metallicity. Some previous stud¬ 
ies have found that Lya escape is easier in low metallic¬ 
ity galaxies, while others have found no dependence on 
metallicity. Gowie et al. (2011) showed that GALEX 
z^O.3 LAEs have lower metallicities than UV-selected 
galaxies with similar UV magnitudes but no Lya emis¬ 
sion. Einkelstein et al. (2011) also showed that GALEX 
z^O.3 LAEs have lower metallicities than similar-mass 
galaxies from the SDSS. In contrast, Atek et al. (2014) 
found that there is no correlation between and 

metallicity for a combined z=0 - 0.3 LAE sample. Hayes 
et al. (2014) showed that in the LARS sample of 
14 LAEs, a few galaxies with high have on av¬ 

erage lower metallicity than the other galaxies with 
low Note that the metallicities in those studies 

are measured from the strong line indexes ([NH]/Ha, 
([OH]-|-[OHI])/H/3, or ([OHI]/H/3)/([NH]/Ha) and have 
large uncertainties. Eor the present Green Pea sam¬ 
ple, we show the relation between and metal- 

licity measured from the more accurate 7). method 
(right panel of figure 6). If both outliers GP0303-0759 
and GP1249-I-1234 are excluded, anti-correlates 
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Fig. 4.— Relations between and velocity quantities of Lya profile. (Upper-left): and velocity of blue peak of Lya profile. 

(Upper-right); and velocity of blue peak normalized by FWHM(Ha). (Lower-left): and velocity of red peak of Lyo profile. 

(Lower-right): and peak separation of Lya profile. The Spearman correlation coefficient r, and null probability P, are given in each 

panel. 


with metallicity with Spearman coefficient r=-0.82 and 
P=0.004 (note that r=-0.58 and P=0.06 if only the AGN 
candidate GP0303-0759 is excluded). It is tighter than 
the correlation between and dust extinction. 

6. Lya RADIATIVE TRANSFER MODELING 

To extract more information from the Lya profiles and 
explore the physical process of Lya escape, we fit the 
Lya profiles of Green Peas with expanding HI shell ra¬ 
diative transfer models (Dijkstra et al. 2014; Gronke et 
al. 2015). In the radiative transfer model, Lya pho¬ 
tons were generated in the shell center and then scat¬ 
tered/absorbed by a shell of HI gas. The intrinsic Lya 
line has a Gaussian profile with width a. The shell is 


described by four parameters: (a) outflow velocity Vexp, 
(b) hydrogen column density Nhi, ( c ) temperature T 
(including turbulent motion as well as the true tempera¬ 
ture), and (d) dust optical depth r^. The radiative trans¬ 
fer simulation generates a grid of Lya spectra covering 
large ranges of parameters. To compare the observed 
spectrum with models, we calculate the and likeli¬ 
hood values, and find the best matched model param¬ 
eters (cr, Vexp, Nhi, T, t^) using Markov Chain Monte 
Carlo (MCMC) methods. The one-dimensional and two 
dimensional marginal posterior distributions are calcu¬ 
lated from the MCMC chains. The resulting param¬ 
eter value and lower and upper uncertainties are the 
50%, 16%, and 84% percentiles respectively in the one- 
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Fig. 5.— Relations between feJc^ and: (Left): the flux ratio of the blue to the red part of the Lyo: profile, 
ratio of the valley to the red peak of the Lyo profile. (Right): the FWHM of the red part of the Lyo profile, 
coefficient r, and null probability P, are given in each panel. 


(Middle): the flux density 
The Spearman correlation 




12+log(0/H) (Tg method) 


Fig. 6. — Relations between feJc^ and: (Left): dust extinction E(B-V). (Right): metallicity. The Spearman correlation coefficient r, and 
null probability P, are given in each panel. GP0303-0759 (marked by a large circle) has the lowest metallicity and E(B-V). Its WISE three 
band colors suggest it is probably an AGN (Malhotra et al. in preparation). GP1249+1234 (marked by a large square) has single peak 
profile and large E(B-V). We exclude both sources when calculating correlation coefficients. 
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Fig. 7. — Comparison of the observed Lyo profiles (blue lines) and the best fit Lya profiles (green lines) from radiative transfer models. 
To show the uncertainty in the fitted profiles, we also plot (grey lines) model Lya profiles corresponding to 25 random samples from the 
MCMC chain of each galaxy. The flux of all profiles are normalized. 


dimensional marginal posterior distribution of each pa¬ 
rameter. The details of the models and fitting methods 
are described in Gronke et al. (2015). As an example, 
we show the marginal posterior distributions and fitting 
result for GP0911-I-1831 in figure 8. 

The model gives good fits for profiles of 9/12 Green 
Peas. For the 9 cases where the model fitting pro¬ 
ceeded smoothly, the best-fit model profiles are shown 
in figure 7, and the derived parameters and uncertain¬ 
ties are shown in Table 4. The best-fit expanding ve¬ 
locities are between 27-127 km s“^ for 8 Green Peas 
except GP1249-I-1234 which has Vexp=354 km s“^ and 
single peak Lya profile. The derived N^j are between 
1019 _ io20cm-2. 

The model fails to fit the three double-peaked profiles 
(GP1424-h4217, GP1133-k6514, GP1219-kl526) where 
the is high and V(Valley)>0 (see Section 4 and 

figure 2). This could be due to a few reasons: (1) ill- 
understood COS instrumental effects in taking spectra 
for slightly extended objects (as the Green Peas); (2) 
inflowing gas; (3) Spatial and velocity offsets of star¬ 
forming regions relative to the Lya emitting gas; (4) 
general inapplicability of simple shell models in some 
cases. We will report fitting with improved model in 
a following paper. For these three cases, we manu¬ 
ally adjust the model parameters to match the observed 
depth of the “valley” and the relative heights of the blue 


and red peaks. The resulting model parameters (Table 
4) show outflowing velocities ^ 20—110 km s“^ and 
JVp[i ~ 10^® ° cm~^. The resulting N^i is moderately 
reliable because a significant larger Nhi can’t fit the ob¬ 
served large residual flux near the line center. 

We plot as a function of fitted Nhi in figure 9. 
For the three cases that needed manual adjusting, we set 
generous errors of Nhi to 0.5 dex. The results suggest an 
anti-correlation between and Nhi with Spearman 

r=-0.69 and P=0.01. Previous studies have suggested 
that LAEs have lower Nhi than non-LAEs (e.g. Shibuya 
et al. 2014; Erb et al. 2014; Hashimoto et al. 2015). Here 
for the first time, we show that is higher at lower 

Nhi in an LAE sample. Therefore we conclude that the 
low column density of HI gas helps Lya escape in these 
Green Peas. 

7. DISCUSSION 

7.1. Correlations between and Lya Profiles 

High quality Lya profiles of Green Peas provide a good 
opportunity to study Lya escape in LAEs. We show that 
the is correlated with kinematic features of Lya 

escape. What are the implications of these correlations 
about Lya escape? 

(1) and V(blue-peak): As shown in Section 4, 
the correlates strongly with V(blue-peak) and peak 
separation. Both V(red-peak) and V(blue-peak) have a 
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Fig. 8.— The fitting results for GP091H-1831. Upper right: The observed Lyo: profile and the best fit model profile. Main plot: The 
one-dimensional and two dimensional marginal posterior distributions from the MCMC chains. In the plot of one-dimensional marginal 
posterior distribution, the dashed lines mark the 16%, 50%, and 84% percentiles (numbers are shown in top of each column). In the plot 
of two dimensional marginal posterior distributions, the blue lines mark the (2, 1.5, 1, 0.5) cr contours, and the gray shading gives the 
posterior density (Gronke et al. 2015). 


complex dependence on gas outflows, and gas tem¬ 
perature. The velocity offset V(blue-peak) is generally 
larger when there is more HI gas and a larger numbers of 
scatterings for Lya photons (e.g. Verhamme et al. 2015). 
So the correlation between and V(blue-peak) sug¬ 
gests LAEs with higher have less scatterings for 

Lya photons and lower Nhi- 
(2) fesT and /a( valley)// a( red-peak): A few 
Green Peas show significant residual flux near Lya line 
center. As the optical depth of HI gas at the systemic 
redshift increases rapidly with Nhi, only the Green Peas 
with low Nhi will have residual flux at Lya line cen¬ 
ter. Low Nhi should correspond to high escape frac¬ 
tion, and indeed we see that is correlated with 

/a (valley) //a (red-peak). 


(3) and FWHM(red): is anti-correlated 

with FWHM(red). To broaden the red peak profile gen¬ 
erally requires scattering of Lya photons to redder fre¬ 
quency by HI gas. Again, more scattering results in more 
dust extinction and smaller 

7.2. Difficulties of Expanding HI Shell Approximations 

We fit the Lya profiles with an expanding single shell 
radiative transfer model. While the models often yielded 
remarkably good fits to the observed line profiles, they 
have limitations that we discuss here. 

Velocity of central valley of Lya profile: The 
single shell radiative transfer model can’t fit the three 
profiles with V{valley) > 0. Simply allowing a negative 
shell velocity (corresponding to a single infalling shell) 
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TABLE 4 

Lyq profile Model Parameters 


ID 

a 

Vecrp 

\og{NHi) 

log(T) 



(km s“^) 

(km s“^) 

(cm~^) 

(K) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

GP1457-I-2232 

280/34 

-12 

20.09/°;22 

4.5/°;l 

'^•'^^-0.28 

GP0303-0759 

245/1^ 

113/11 

19-4l/°o;(,l 

'^•'-0.5 

u.zo_Q_i8 

GP1244-I-0216 

351/1/ 

55 //^ 

19.80/°;°/ 

4 O+0.3 

o-i4/S;ll 

GP1054-I-5238 

226/)/ 

91/1° 

19.78/°;°/ 

4 c+0.2 

0.54/°;f 

GP1137-I-3524 

239/(2 

127/^3 

19.60/°;°? 

0 ^-1-0.2 

Q 47-^0.14 

GP0911-I-1831 

215 /® 

120/1 

19.00/°;°/ 

0 ^-1-0.2 

f) kq+0.12 

D.OU-0.09 

GP0926-I-4428 

266/1 

110/1 

19.40/°;°? 

q 7'-|-0.2 
' -0.6 

D.UU-0.00 

GP0815-I-2156 

207/«2 

34 /I 3 

19.23/°;// 

4 q+0-2 

^•^-0.3 

f) ni +0.25 

GP1249-1-1234 

206/^“ 

354//1 

19 fil 

iy.Di o Qg 

^•^-1.0 

-*^•^^-0.14 

GPT424+42T7 

350 

no 

18.0 

3.4 

0.16 

GP1133-I-6514 

214 

20 

17.8 

3.0 

0.16 

GP1219-I-1526 

214 

20 

18.0 

3.0 

0.16 


Note. — Column Descriptions: (2) Icr width of the Gaussian profile of intrinsic Lyo; line; (3) outflowing velocity of HI gas shell; (4) HI column 
density of the expanding gas shell; (5) temperature of HI gas; (6) dust optical depth. The first 9 rows show the best fit parameters and errors for 
nine Green Peas. The 68% uncertainties are given for all parameters. The last 3 rows show the manually adjusting results for three profiles that 
the model fails to fit. 



(cm 2 ) 


Fig. 9.— Relations between /i^c^ and best fit Nj{j from the 
expanding shell model. N}ji was derived by fitting of the shell 
model parameters for 9/12 Green Peas. For the other three cases 
(GP1424+4217, GP1133+6514, and GP1219+1526, marked by 
large circles) where the fitting procedure failed, we plot the Njji 
obtained by manually adjusting the model parameters to match 
the observed depth of the “valley” and the relative heights of blue 
and red peaks (see Section 6). GP1249+1234 (marked by a large 
square) has single peak profile, large expanding velocity, and large 
E(B-V). The Spearman correlation coefficient r, and null probabil¬ 
ity P, are calculated with all 12 sources. 

would permit V{valley) > 0, but is not a solution, be¬ 
cause it also produces profiles whose blue peak is higher 
than the red peak (contrary to observation). One possi¬ 
bility is that besides the expanding HI shell component, 
there is another inflowing HI component that absorbs 
more Lya photons in the red side of profile and make 
V{valley) > 0. Another possibility is that there are a 
few star-forming gas clumps with different central veloc¬ 
ities. Thus there could be spatial and velocity offsets of 
star-forming regions relative to the Lyo emitting gas. If 
the Lya emissions are dominated by the gas clumps with 
redshifted velocities relative to the average systemic red- 


shift, then the observed total Lya profile would show 
a red-ward offset and V{valley) > 0. Generally, to fit 
these profiles requires us to consider more realistic HI 
gas distributions than the HI shell approximation. 

Lya intrinsic line width; The model fitting results 
show intrinsic Lya line width between 200 - 350 km s“^ 
(1 a of Gaussian profile) and FWHM about 470-820 km 
s“^, which are ^ 2-3 times larger than the intrinsic Ha 
line width. If the intrinsic Lya profile is as narrow as 
Ha profile, producing the observed broad Lya profiles 
would require large N^i and a high HI gas temperature 
T. However larger N^j and T will also increase absorp¬ 
tion of Lya photons, and result in strong central absorp¬ 
tion and weak blue peak, which contradict the observed 
profile. Hashimoto et al. (2015) fit Lya profiles of a few 
z^2 LAEs and also found the best fit intrinsic Lya width 
is a few times larger than the intrinsic Ha or [OIHJA5007 
width. Lya and Ha powered by star formation ought to 
originate in the same gas, so a wider intrinsic Lya profile 
would suggest either that not all the Lya is powered by 
star formation (cf. Hashimoto et al. 2015), or that there 
are important radiative transfer effects that broaden Lya 
nearer to the source, before it encounters the model’s ex¬ 
panding outer HI shell. 

Shell expansion velocity: The resulting shell ex¬ 
pansion velocities are between 20-130 km s“^ for 11/12 
Green Peas (and larger for the single-peaked line in 
GP1249-I-1234). These expansion velocities are compara¬ 
ble to the 1(7 width of Ha emission lines (80-130 km s“^), 
but are about 2 times smaller than the average outflow 
velocities traced by low ionization metal absorption lines 
(^ 80-300 km s“ , see Henry et al. 2015). This may 
suggest that the low ionization metal absorption lines 
trace a different outflowing gas component from the HI 
gas, and that the velocity distribution of HI gas is more 
similar to that of HII regions. 

7.3. Metallicity and Lya Escape 

We showed that anti-correlates with metallicity. 
The /iir is also higher in Green Peas with lower stellar 
mass and higher sSFR. However, only dust and HI gas 
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directly interact with Lya photons. What causes the cor¬ 
relation of and metallicity? It is not simply caused 
by a correlation between the abundances of metal and 
dust, as the metallicity has tighter correlation with 
than dust reddening. Thus the metallicity corre¬ 

lation may suggest lower Nhi in lower metallicity Green 
Peas. As the mass and metallicity are correlated, one 
possibility is galaxies with lower mass and lower metal¬ 
licity have less HI gas. The other possibility is outflows 
in lower metallicity galaxies blow out more holes with 
low Nhi and increase 

Outflows have a large impact on the metallicity of low 
mass galaxies. These emission line selected galaxies (in¬ 
cluding LAEs and Green Peas) can have metallicities as 
low as 1 dex below the usual mass - metallicity rela¬ 
tion (e.g. Xia et al. 2012; Ly et al. 2014; Song et al. 
2014). The galactic mass - metallicity relation is gener¬ 
ally a result of interactions between the metal enrichment 
and gas outflows: star formation enriches the interstel¬ 
lar medium, but outflows can drive metal enriched gas 
out and decrease the metallicity. Since outflows are more 
effective where the potential is shallow, this process re¬ 
sults in lower metallicities at lower masses. Emission line 
galaxies like LAEs and Green Peas tend to have compact 
star formation and high sSFR, resulting in stronger out¬ 
flows and lower metallicity than other galaxies with sim¬ 
ilar mass. The outflows from stellar feedback are dom¬ 
inated by ionized gas. These ionized outflows can drive 
HI gas out or make ionized holes in interstellar medium, 
thus reducing the effective column of Nhi for Lya pho¬ 
tons. 

8 . CONCLUSION 

We have investigated high quality Lya spectra of a 
sample of 12 Green Peas. All show Lya emission lines. 
We compared the EW(Lya) distribution of this Green 
Peas sample to high redshift LAEs. From the Ha emis¬ 
sion lines we obtained the systemic redshift and intrinsic 
Lya emission. The high S/N Lya spectra permit us to 
measure accurately the blue peaks of Lya. So we mea¬ 
sured the velocities and width of the Lya profiles and 
explored correlations of profile features with Lya escape 
fraction. We also compared with dust reddening 

and metallicity. In addition, we fit Lya prohles with ra¬ 
diative transfer model and discussed constraints on HI 
gas and Lya escape. Our main results are as follows. 

1. The EW(Lya) distribution of this Green Pea sam¬ 
ple is very similar to our z=2.8 LAE sample. As 
these Green Peas have small mass, compact size. 
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and high emission line equivalent widths, similar 
to high-z LAEs, we conclude that Green Peas are 
the best analogs of high-z LAEs in the local uni¬ 
verse. 

2. The shows correlations with the blue peak ve¬ 
locity of Lya, the ratio V(blue-peak)/FWHM(Ha), 
the flux density ratio /a( valley)//^(red-peak), and 
the width of Lya red peak FWHM(red). As more 
scatterings in HI gas can make the Lya blue peak 
bluer, the residual central flux smaller, and the red 
peak profile wider, these correlations strongly sug¬ 
gest low Nhi and fewer scatterings help the Lya 
photons escape. 

3. The single shell radiative transfer models can re¬ 
produce most profiles, and get column densities 
of HI gas. The resulting Nhi anti-correlates with 

in this LAE sample, again indicating that low 
Nhi is key to Lya escape. 

4. However, the single shell model also has difficul¬ 
ties. It fails to ht three profiles. The best-ht in¬ 
trinsic line width is 2-3 times larger than the Ha 
width. The best-fit HI gas outflow velocity is small 
and only about half of the average velocity of low 
ionization absorption lines. Fully reproducing all 
features of the observed lines will require model¬ 
ing Lya radiative transfer in more realistic HI gas 
distributions. 

5. In this Green Pea sample, shows a weak anti¬ 
correlation with dust reddening, and a stronger 
anti-correlation with metallicity. We suggest that 
this correlation may be an effect of the mass- 
metallicity relation and trends of other properties 
with galaxy mass. Lower metallicity galaxies are 
likely to have less HI gas and dust. In addition, ion¬ 
ized gas outflows can blow out the metal enriched 
gas, making holes with low Nhi and helping Lya 
escape. 

In conclusion, Green Peas provide an unmatched op¬ 
portunity to study Lya escape in LAEs. Our results 
suggest that LAEs with high Lya escape fraction have 
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outflow. In future work, we will compare Lya profiles 
and the correlations found in this work with improved 
models and provide more quantitative constraints on the 
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